Stability of charged solitons and formation of boson stars in 
5-dimensional Anti-de Sitter space-time 

Yves Brihaye **, Betti Hartmann ' and Sardor Tojiev * 



^Physique-Mathematique, Universite de Mons-Hainaut, 7000 Mons, Belgium 



CO ■ 1 School of En gineering and Science, Jacobs University Bremen, 28759 Bremen, Germany 

O : January 14, 2013 

a 

a 

Abstract 

We study the stability of charged solitons in 5-dimensional Anti-de Sitter (AdS) space-time. We show 
that for appropriate choices of the parameters of the model these solutions become unstable to form scalar 
hair. We find that the existence of charged solitons with scalar hair depends crucially on the charge and 
the mass of the scalar field. We investigate the dependence of the spectrum of solutions on the mass of the 
scalar field in detail. For positive mass of the scalar field the hairy solitons can be interpreted as charged 
boson stars. We find that for sufficiently small value of the charge of the scalar field a "forbidden band" of 
the boson star mass and charge exists, while all our results indicate that - contrary to the asymptotically 
flat space-time case - boson stars in asymptotically AdS can have arbitrarily large charge and mass. 
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t— I , The Anti-de Sitter (AdS)/Conformal Field theory (CFT) correspondence relates gravity theories in (d + 1)- 
' dimensional asymptotically AdS space-time to a CFT living on the d-dimensional boundary of that space-time 
, [TJ [2]. As such, classical solutions in asymptotically AdS have gained a lot of interest. This includes both 
. . ' black hole solutions as well as globally regular, solitonic-like solutions. Within the context of the holographic 
^ , description of high-temperature superconductivity solutions with asymptotic planar AdS have been considered 
' [31 IH [H |5] . It was shown that the formation of scalar (or vector) "hair" on the solutions is the dual description 
5_i , of the onset of superconductivity. This is possible since the effective mass of the scalar field drops below the 
. . .' Breitenlohner-Freedman (BF) bound [7] under certain circumstances and the black holes or solitons, respectively, 
become unstable to the formation of scalar hair possessing however still an asymptotic AdS space-time. Besides 
being interesting due to the holographic interpretation, the stability of classical field theory solutions is, of course, 
also of interest by itself as it could shed well light on other interesting questions, e.g. black hole uniqueness. 
Consequently, the stability of black holes and solitons in global AdS as well as in asymptotic hyperbolic AdS 
have also been discussed. 

In [8] uncharged black holes in (4 + 1) dimensions have been considered. It was shown that static black holes 
with hyperbolic horizons can become unstable to the formation of uncharged scalar hair on the horizon of the 
black hole due to the existence of an extremal limit with near- horizon geometry AdS 2 x H 3 0(101 E]. These 
studies were extended to include higher order curvature corrections in the form of Gauss-Bonnet terms [12j . 
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Charged black hole and soliton solutions in asymptotically global AdS in (3 + 1) dimensions were studied 
in [13] . It was pointed out that the solutions tend to their planar counterparts for large charges since in 
that case the solutions can become comparable in size to the AdS radius. The influence of the Gauss-Bonnet 
corrections on the instability of these solutions has been discussed in [14]. The corresponding investigation in 
(4 + l)-dimensional global AdS was done in [HI [16]. The existence of solitons in this case had been suggested 
previously in a perturbative approach [17]. In all studies the mass of the scalar field had been set to zero. In 
[18] these results were extended to a tachyonic scalar field as well as to the rotating case. Recently, charged 
soliton solutions with positive scalar field mass have been studied in (3 + l)-dimensional global AdS [T5]. In 
this case, the solutions carrying scalar hair can be interpreted as charged non-spinning boson star solutions in 
global AdS space-time. Uncharged boson stars in AdS were first discussed in (d + l)-dimensional AdS space- 
time using a massive scalar field without self-interaction |20] and with an exponential self-interaction potential 
[2"Tl |2"2"] . respectively. Spinning solutions in (2 + 1) and (3 + 1) dimensions have been constructed in [S3] and 
[2~4") . respectively. 

In this paper, we arc interested in the stability of charged solitons in (4 + l)-dimcnsional global AdS with 
respect to formation of scalar hair. We investigate in detail the dependence of the pattern of solutions on the 
scalar field mass and charge considering negative and positive values for the mass alike. The work extends the 
analysis of [15] and of [25] to now include positive values for the mass such that the solutions with scalar hair 
in our case can be interpreted as charged boson stars in (4 + l)-dimensional global AdS. 

This paper is organized as follow: in Section 2, we present the model, give the ansatz, the equations of 
motion and the relevant boundary conditions. It is well known that the generic (hairy) solutions cannot be 
constructed in terms of elementary functions and that numerical techniques are needed. The discussion of the 
numerical solutions is given in Section 3. Finally, Section 4 contains the conclusions. 



2 The model 

In this paper, we are studying the formation of scalar hair on electrically charged solitons in (4 + l)-dimensional 
Anti-de Sitter space-time. The action reads : 

S = J d 5 x^ {R-2A + 16^G£ mattcr ) , (1) 

where A = — 6/L 2 is the cosmological constant, R denotes the Ricci scalar and G is Newton's constant. £ ma tter 
denotes the matter Lagrangian which reads : 

£ matter = --^F MN F MN -{D M ^Y D M ^-m 2 r^P , M,N = 0, 1,2,3,4 , (2) 

where Fmn = QmAn — OnAm is the field strength tensor and Z?m"0 = c'mV' - ieAj^tp is the covariant derivative, 
e and m 2 denote the electric charge and mass of the scalar field ip, respectively. 
We choose the following spherically symmetric Ansatz for the metric : 

ds 2 = -f(r)a 2 (r)dt 2 + -±-dr 2 + ^dtt 2 , (3) 
f(r) L 2 

where / and a are functions of r only and dfl 2 denotes the line element of a unit 3-dimensional sphere. 
For the electromagnetic field and the scalar field we have [4] : 

A M dx M = (j}{r)dt , ^ = VKO (4) 

such that the solutions possess only electric charge. 

The coupled gravity and matter field equations are obtained from the variation of the action with respect 
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to the matter and metric fields, respectively, and read 



/' = ^l-/ + 2^)- 7 ^(2eW + /(2m 2 aV + 0' 2 ) + 2/ 2 aV 2 ) , (5) 
r(eW+q 2 /y 2 ) 

a = 7 ' ( 6 ) 

3 a '^"2^, (7) 



raj f 

-( 3 Z + f 4 + ^)r'-[^--)r. us, 



r / a/ V/ 2fl2 / 

where 7 = 167rG. Here and in the following the prime denotes the derivative with respect to r. These equations 
depend on the following independent constants: Newton's constant G, the cosmological constant A (or Anti-de 
Sitter radius L) and the charge e and mass m of the scalar field. 
The system possesses two scaling symmetries: 

r -> Ar , t-t Xt , L — s- XL , e — > e/A (9) 

as well as 

-> A0 , V ^ , e e/A , , (10) 

which we can use to set L = 1 and 7 to some fixed value without loosing generality. 

In order to solve the equations, we have to fix appropriate boundary conditions. For soliton solutions regular 
on r G [0, 00 [ we require 

/(0) = 1 , 0'(O) = O , V(0) = ^o , V>'(o) = o , (11) 

where "00 is a f ree parameter. 

Asymptotically, we want the space-time to be that of global AdS, i.e. we can choose a(r —> 00) — > 1. Other 
choices of the asymptotic value of a(r) would simply correspond to a rescaling of the time coordinate. The 
matter fields on the other hand obey : 

*(r>l) = /*-£ , ^r»l) = ^ + ^ (12) 

with 

A_ = 2 - ^/A + m 2 L 2 , A+ = 2 + ^4 + m 2 L 2 . (13) 

The parameters (1, Q are the chemical potential and the electric charge, respectively. Note that for m 2 > we 
have to choose -0- = 0. We will hence also adapt this choice for m 2 < 0. Within the holographic interpretation 
of our model ip+ will correspond to the expectation value (O) of the operator O dual to the scalar field on the 
conformal boundary of AdS. Note that in comparison to the case of holographic superconductors (see e.g. [4] 
and references therein) the space-time in this study possesses spherically symmetric sections for constant r and 
t. In other words, for r — > 00 our space-time corresponds to global AdS. 

The mass of the solution is given in terms of the coefficients that appear in the asymptotic fall-off of the 
metric functions which reads 

/(r»l) = l + ^ + ^ + ... , o(r>l) = l + £ + ... . (14) 

fi and C4 are constants that have to be determined numerically and that depend on the couplings in the model. 
The mass M then reads 

M = -y- ( 15 ) 
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3 Numerical results 



We have used a Newton- Raphson algorithm with adaptive grid scheme [26] to construct the solutions numerically. 
In the numerical construction, we set L 2 = 1 and 7 = 9/40 without loss of generality. The remaining parameters 
of the model are e 2 and ra 2 . 

Once fixing e 2 , m 2 families of regular solutions labelled by ip(0) = 0o (and then by Q) can be constructed 
numerically. The pattern of solutions turns out to be so involved that none of these quantities can be used to 
parametrize the branches uniquely: in some regions two or more solutions are characterized by the same value 
of Q or of ipo . Furthermore, for the values of the parameters studied we noticed the occurrence of at least two 
disconnected branches. A similar phenomenon was observed in different models previously [151 127) . 

3.1 m 2 = 1 

We first discuss the solutions for positive values of the mass parameter m 2 choosing m 2 = 1. In this case the 
soliton solutions carrying scalar hair can be interpreted as charged boson stars in global AdS. 

Fixing in addition e 2 to a particular value we have constructed a first branch of solutions which can be 
labelled by ipo = i/j(0). Letting -00 -> we find that 0(r) = ip(r) = 0, while a(r) = 1. The resulting solution 
hence corresponds to global AdS space-time. Accordingly, the mass M and the charge Q tend to zero in this 
limit. 

We find that the solutions' properties are crucially related to the choice of e 2 . We first choose three specific 
value of e 2 to demonstrate the general pattern. Our results are shown in FigfTJ where we give M and a(0) 
as functions of '0(0) = 0o (upper and lower figure, respectively) and also present M as function of Q (middle 
figure). Note that the middle and lower figure do not contain the results for e 2 = 3. This is due to the fact that 
the data for e 2 = 3 is so close to that for e 2 = 2 that the curves cannot be distinguished. Furthermore, in the 
M-over-Q plot the curve for e 2 = 3 would be continuous and hence the gap for e 2 = 2 and e 2 = 1 would not 
be clearly visible. The solid lines correspond to the fundamental solutions for which the scalar field function 
ip(r) decreases monotonically from 0o > to zero at r = 00. On the other hand the dashed line (for e 2 = 1) 
corresponds to first excited solutions, i.e. solutions that possess one node in the scalar field function. 

Our results suggest the existence of a critical value e 2 r of the parameter e 2 that separates the pattern of 
solutions into (at least) two types. To be more precise, we find that for e 2 < e 2 r two branches of solitons with 
scalar hair can be constructed. One branch - which we will call the "main branch" in the following - is connected 
to the AdS vacuum at "0(0) = (and M = Q = 0, <z(0) = 1) and can be continued to arbitrarily high values of 
ipo (see label A in FigfT]). In this limit, the mass M and Q tend to finite values, while a(0) tends to zero. Note 
that in this case the so-called "planar limit" [37] which appears for large Q cannot be reached by a smooth 
deformation of the AdS vacuum. Next to this main branch, a second branch of hairy soliton exists (label B and 
label C in FigflJ. This is completely disjoint for the main branch. Again, solutions exist for arbitrarily high 
values of 0o- In this limit (see label B) the mass M and the charge Q are again finite, however higher in value 
than the corresponding values in the same limit on the main branch (see label A). All our numerical results 
indicate that the two branches do not join for sufficiently large values of ipo, but remain more or less parallel. 
This leads to the observation that a "forbidden band" for the mass M and equivalently for the charge Q exists 
for charged boson stars in asymptotically AdS if the charge of the scalar field e 2 is "small enough" . 

In the limit 0o - > the solutions on the second branch do not tend to the AdS vacuum, but the mass M 
and charge Q tend to infinity (see label C). All these features are clearly visible in FigfTJfor e 2 = 2. For e 2 = 1 
we also find the main branch and all our numerical results indicate that only this branch exists. If the second 
branch also existed it would be extremely difficult to construct this numerically. 

Another striking feature of the second branch of solutions (see label C in Fig. []} is that it apparently 
extends backwards to ipo = although the numerical construction becomes difficult in this limit. A solution 
in this parameter range is given in Fig. [2] We observe that the metric function f(r) possesses a minimum at 
some intermediate value of the radial coordinate r = r m ; n which increases with decreasing ipQ. Furthermore, 
the scalar field ip(r) deviates only little from the value tpa up to r = r m i n , while for r > r m - m the scalar function 
reaches its power-like behaviour. 

Let us also remark that the qualitative features are similar for the first excited solutions, i.e. the solutions 
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Figure 1: We show the mass M (upper figure) and the value of a(0) (lower figure) as functions of ip$ for m 2 = 1 
and for e 2 = 1 (black), e 2 = 2 (red) and e 2 = 3 (blue), respectively. We also give the mass M in dependence 
on Q (middle figure). Note that the middle and lower figure do not contain the results for e 2 = 3. This is 
due to the fact that the data for e 2 = 3 is so close to that for e 2 = 2 that the curves cannot be distinguished. 
Furthermore, in the M-over-Q plot the curve for e 2 = 3 would be continuous and hence the gap for e 2 = 2 and 
e 2 = 1 would not be clearly visible. In all three figures the solid lines represent the fundamental solutions, while 
the dashed lines correspond to the first excited solutions (here given for e 2 = 1). 
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Figure 2: We show the profiles of the metric functions f(r) (solid) and a(r) (long-dashed) as well as of the 
scalar field function ijj(r) (short-dashed) for e 2 = 2, m 2 = 1 and for ip(0) = 1.5 (black) and ip(0) = 0.9 (red), 
respectively. 

that possess one zero in the scalar field function. This is shown for e 2 = 1 in Fig Q] (dashed line). The dependence 
of the mass M on -0o is similar than that of the fundamental solution and the M-Q-plot also clearly shows that 
the mass of the excited solutions is larger than that of the fundamental solution for a fixed value of Q. These 
solutions are hence unstable to decay to the fundamental ones. 

For e 2 > e 2 r we find that only one branch of solutions exists. This branch can be constructed from V'(O) = 
(and M = Q = 0) at the AdS vacuum up to a maximal finite value of ip(0) = '(/'(0) max . From this maximal 
value of V'(O) the branch can be continued when decreasing ^(O) again. On this back-bending branch the mass 
M and charge Q of the solution tends to infinity. Accordingly, the solutions tend to the "planar limit" [37] in 
which the size of the solitons becomes comparable to the AdS radius. The shape of the solution on this branch 
is reminiscent of the so-called "thick-wall limit" [28 of uncharged Q-balls and boson stars which appears for 
ipo — > 0. This also means that we can find a branch of solutions that interpolates smoothly between the AdS 
vacuum and the planar limit. Finally, we find that for e 2 = e 2 r the main branch and the second branch join at 
a specific value of ip(Q). 

We have then done a detailed analysis of the solutions around the critical value of e 2 . Our results are given 
in Fig|3l where we show the mass M as function of ip(0) for several values of e 2 close to the critical value of e 2 
which we determine to be e 2 r sa 2.6529 for m 2 = 1. This should be compared to the value for m 2 = 0, which is 
e 2 r sa 2.3696. Hence, e 2 r increases with increasing m 2 . 

These results can be compared to the analogue case in (3 + 1) dimensions [T§] . 

3.2 m 2 = -3 

If we want the asymptotic AdS of our space-time to be stable the mass to 2 of the scalar field is limited by the 
BF bound (with our choice L = 1 we have m gF — —4 in (4 + 1) dimensions). We have therefore studied the 
spectrum of the hairy soliton solutions for a value of m 2 relatively close to the BF bound, namely m 2 = —3. 
To appreciate the difference with respect to the case m 2 = 1 we present the counterpart of FigJT] in FigJH 
Here, we find that there is again a critical value of e 2 that separates the pattern into distinct types. Our 
numerical results indicate that e 2 r sa 1.3575. For e 2 < e 2 r (here e 2 = 1, black curves) we again find a main 
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Figure 3: The mass M is shown as functions of ^(0) = -00 for to 2 = 1 and several values of e 2 close to the 
critical value e 2 r . 

branch of solutions which is connected to the AdS vacuum at ip(Q) = 0, M — Q = 0. On this main branch 
the solutions can be constructed for arbitrarily large values of ip(0) and in the limit f/'(0) — > oo (label A in 
FigHJ) the mass M and charge Q tend to finite values, while a(0) tends to zero. Similarly, a second branch of 
solutions exists in this case, however, this is only present for ip(0) larger than a minimal value V'min(O) (label 

B) . At this minimal value a third branch of solutions can be constructed that extends back to i/j(0) — > oo (label 

C) . We hence find that for ip(0) > V'min(O) up to three solutions with different masses exist for the same value 
of tp(0). However, when considering the physically relevant parameters M and Q the solutions are uniquely 
characterized by these values (see M-Q-plot in Fig0]). Again, we find a "forbidden band" of the mass M and 
charge Q (between label A and label B). 

For e 2 > e 2 r (here e 2 = 2, red lines) we find that only one main branch of fundamental solutions can be 
constructed. Contrary to the case with positive m 2 this exists for all values of ipo > with the mass M diverging 
for ipo — > oo. In this case, solutions exist for all possible values of the mass M and charge Q. 

We were also able to construct a branch of excited solutions for e 2 > e 2 r for which the scalar field function 
possesses a node. This is indicated by A' in FigHJ These solutions exist only down to a minimal value of the 
mass M and charge Q. 

3.3 The general pattern 

After choosing particular values of the scalar field mass m 2 we would now like to study the general pattern. For 
all values of e 2 and fixed, but arbitrary m 2 we find that a main branch of solution exists that is connected to 
global AdS in the limit ipo 0. For e 2 < e 2 r , the charge Q and mass M of the solutions on this branch tend 
to finite values for ip — ¥ oo and a second, disjoint branch of solutions exists. This leads to the observation that 
in these cases a "forbidden band" of the mass M and charge Q exists. On the other hand for e 2 > e 2 r the mass 
M and charge Q can have arbitrarily large values in this limit. We have observed that e 2 r depends on m 2 and 
increases with to 2 . Our numerical results indicate the approximate relation 

e 2 r «2.4+^. (16) 
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Figure 4: We show the mass M (upper figure) and the value of a(0) (lower figure) as functions of V>o for m 2 = —3 
and for e 2 — 1 (black) and e 2 — 2 (red), respectively. We also give the mass M in dependence on Q (middle 
figure). In all three figures the labels A, B, C represent branches of fundamental solutions, while A' corresponds 
to a branch of first excited solutions (here given for e 2 = 1). 
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3.4 Dependence on m 2 

As we have seen above, the M-Q-plot looks qualitatively similar for all possible values of m 2 , however the 
properties of the solutions given as function of ipo depend on m? . We have hence studied the properties of the 
solutions in dependence on m 2 in detail by fixing and e 2 and varying m 2 € [—3, 1]. 

We find that for those values of ipo an d e 2 for which only one solution exists (which is then a solution on 
the branch connected to the AdS vacuum) a variation of m 2 leads to a smooth deformation of the solutions. 
However, we have found it numerically more challenging to vary m 2 for those solutions that are on branches 
not connected to the AdS vacuum. The reason is that in this case several solutions for one fixed value of "00 
co-exist. 

Let us discuss one example of this type in the following. We have chosen branch A' given in Fig. |U This 
is a branch of first excited solutions possessing one node in the scalar field function. Moreover, the metric 
function f(r) possesses a local minimum at say r = r m - ln . We find that fixing ipQ and increasing m 2 the value 
of f(r = r min ) decreases and tends to zero at some critical value of m 2 = m 2 T . This is shown in FigjSJ where 
we give the metric functions (upper figure) and matter functions (lower figure) for e 2 = 2, ?/>o = 5-0 and two 
different values of m 2 . In this case, we find f(r = r m j n ) « for m 2 w 0. At the approach of this critical point 
we find that ip(r > r m j n ) m and a(r > r m ; n ) w 1 suggesting that the limiting solution is an extremal Reissner- 
Nordstrom-AdS solution with horizon at r = r m j n . This is confirmed by the observation that the functions f(r) 
and 4>{r) are given by the analogue expressions for the given values of M and Q. For r < r m ; n the scalar field 
remains non-trivial, while the electric potential <f>{r < r m i n ) w 0. This suggests that the space-time possesses 
two distinct regions matched at r — r m i n : for r < r m j n the space-time is that of a compact uncharged hairy 
soliton which can be interpreted as an uncharged boson star for m 2 > 0, while for r > r min the space-time is 
given by an extremal Reissner-Nordstrom-AdS solution. As a consequence not all of the solutions occurring for 
negative values of m 2 have a counterpart for m 2 positive. 

All the phenomena are summarized in FigJBl where we show the mass M, the charge Q and 4>(0) as a function 
of m 2 for e 2 = 2 and ip(0) = 3.5. We find that for positive m 2 a branch of solutions exists down to a minimal 
value of m 2 where it bends backwards on a second branch of solutions. The solutions on this second branch 
have higher mass M and charge Q as the solutions for the same value of m 2 on the first branch. These are the 
two disjoint branches labelled A and B in FigfTJ Apparently, these two branches join for sufficiently small m 2 
such that for m 2 small enough only one branch persists (see Figjl]). 

Moreover, for negative m 2 close to the BF bound only one branch of solutions exists that extends up to 
a maximal value of m 2 . At this value of m 2 we find that </>(0) = such that the limiting solution is that 
described above (see FigEJ consisting of an external extremal Reissner-Nordstrom-AdS solution and an internal 
uncharged hairy soliton. 

4 Conclusion 

In this paper, we have studied the formation of scalar hair on charged solitons in global AdS space-time. We 
have studied the dependence of the solutions on the choice of the charge e 2 and the mass m 2 of the scalar 
field. For positive m 2 we can interpret these solutions as charged boson stars in AdS. We find that the pattern 
of solutions depends crucially on the choice of e 2 and m 2 with a critical value e 2 r « 2.4 + m 2 /3 dividing this 
pattern into distinct types. Interestingly, we observe that boson stars in AdS can have arbitrarily large mass 
and charge, however also that a "forbidden band" of the mass and charge at intermediate values of the mass 
and charge exists. 

In the context of the AdS/CFT correspondence solitons in AdS have been used to describe insulator- 
superconductor phase transitions within the condensed matter application [29j [30], but also to describe 
confinement-deconfinement phase transitions in Quantum Chromodynamics [3 1 1 132 ] . In particular, in [30] it 
has been suggested that uncharged, planar boson stars in AdS could be the dual of glueball condensates. In our 
case, the boson stars have an asymptotically global AdS and are charged. It would be interesting to understand 
what the holographic interpretation of our solutions corresponds to. This is currently under investigation. 
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Figure 5: We show (upper figure) the metric functions f(r) (black) and a(r) (red) and (lower hgure) the matter 
functions ip(r) (red) and (p(r) (black) for = 5.0, e 2 = 2 and for m 3 = —3 (dashed lines) and m 2 — (solid 
lines), respectively. 
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